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Hard X-ray surveys are an important tool for the study of active galactic nuclei (AGN); they 
provide almost an unbiased view of absorption in the extragalactic population, allow the study 
of spectral features such as reflection and high energy cut-off which would otherwise be unex- 
plored and favour the discovery of some blazars at high redshift. Here, we present the absorption 
properties of a large sample of INTEGRAL detected AGN, including an update on the fraction 
of Compton thick objects. For a sub-sample of 87 sources, which represent a complete set of 
bright AGN, we will discuss the hard X-ray (20-100 keV) spectral properties, also in conjunction 
with Swift/BAT 58 month data, providing information on BAT/IBIS cross-calibration constant, 
average spectral shape and spectral complexity. For this complete sample, we will also present 
broad-band data using soft X-ray observations, in order to explore the complexity of AGN spectra 
both at low and high energies and to highlight the variety of shapes. Future prospects for AGN 
studies with INTEGRAL will also be outlined. 
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1. Introduction 

The study of active galactic nuclei (AGN) above lOkeV is essential if one wants to study non- 
themial processes and observe those sources which are strongly affected by absorption in the soft 
X-ray band. Another advantage is the possibility of having information on spectral features such 
as the high energy cut-off and the reflection fraction, which cannot be explored with observations 
performed below lOkeV. The determination of these parameters is important for many reasons: 
they provide an insight into the physical properties of the region around the central power source, 
play a key role in synthesis models of the cosmic X-ray background and are important ingredients 



for unification theories and torus studies (e.g [B3, 10, 11 1). 



2. The sample 

In the 4* survey [Q] made by INTEGRAL/IBIS [||] there are 234 objects which have been 
identified with AGN. To this set of objects, one can add 38 galaxies listed in the INTEGRAL all- 
sky survey by [ p^ ] updated on the website' but not included in the Bird et al. catalogue due to 
the different sky coverage. The final dataset comprises 272 AGN (last update March 2011), which 



represents the most complete view of the INTEGRAL extragalactic sky to date [21]. This sample 
has two great strengths: a) all sources have optical spectra, which means a secure identification 
and a measured redshift (the only exception is the BL Lac object RX J0137.7-I-5814 for which 
no redshift is available), and b) all sources have X-ray data available, which provides a measure 
of the source intrinsic absorption. In the left panel of figure |l|, the whole sample is reported in 
the classical 20-100 keV luminosity v^^ redshift plot. The luminosities have been calculated for 
all sources assuming Ho=71 kms^' Mpc^' and qo = 0. From this figure, it can be estimated that 
our sensitivity limit is around l.SxlO^^'ergcm^^ s^'. We find that the source redshifts span from 
0.0014 to 3.7 with a mean at z=0.1477, while the Log of 20-100 keV luminosities ranges from 
40.7 to ~48 with a mean at around 10"*^ ergs"^ NGC 4395 (a Seyfert 2) is the closest and least 
luminous AGN seen by INTEGRAL, while IGR J225 17-1-22 18 (a broad line QSO) is the farthest 
and most luminous; the former hosts a relatively small central black hole (M~10'^-10^ M©), while 
the latter houses a more massive one (M=10^ M0). Thus the INTEGRAL AGN catalogue spans 
many order of magnitudes in source luminosity and black hole mass, indicating that it covers a 
large range in Eddington luminosities and accretion rates. 

Although useful for statistical and population studies, this catalogue is too large and contains 
too many weak sources to be appropriate for hard X-ray spectral studies (either individually or 
in combination with low energy observations). For this purpose we have created a subsample of 
objects which constitute a complete sample of AGN selected in the 20-40 keV band. Full details on 



the extraction of this complete sample are given in [20]. It consists of 87 active galaxies of various 
optical classification: 41 type 1 AGN (Seyfert 1, 1.2 and 1.5), 33 type 2 AGN (Seyfert 1.9 and 2), 5 
narrow line Seyfert Is (NLSy Is) and 8 Blazars (QSOs and BL Lacs). In the right panel of Figure |l], 
the objects in this subsample are displayed in the 20-100 keV luminosity vs redshift plot; it is clear 
from a comparison between the two panels that the data sets are similar in parameter distribution 



http://hea.iki.rssi.ru/integral/survey/catalog.php 
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Figure 1: Left Panel: Hard X-ray luminosity vs redshift for the INTEGRAL sample of 272 AGN. Circles 
are type 1 AGN, squares are type 2 AGN and triangles are blazars. Filled symbols represent objects with 
absorption measured while open symbols refer to sources with upper limits Nh- Right Panel: Same as in the 
left panel but using only objects belonging to the complete AGN sample. 

SO that we can safely assume that the complete sample is representative of the entire population of 
AGN selected in the hard X-ray band. 

3. Absorption properties 

In the following, Nh=10^^ cm^^ is assumed as the dividing line between absorbed and unab- 
sorbed sources: this is the value conventionally used because it corresponds to a column density 



sufficiently high to hide an AGN broad line region (BLR, [|29||). We also consider as Compton thick 
an object with a column density in excess of 1.5 x 10^^ cm^^. It is worth noting that for a number 
of objects we did not measure any absorption in excess to the Galactic value and have therefore 
used the Galactic column density as an upper limit to the source intrinsic absorption. With the 
assumptions made above, the fraction of absorbed objects present in our sample is 49%, while that 
of Compton thick objects is around ~6%, in full agreement with previous estimates available in 
the hterature 

Despite the fact that hard X-ray instruments are the least biased in terms of detecting absorbed 
AGN, they still miss some Compton thick objects, essentially those with weak (intrinsic) fluxes and 
at large distances. This has been fully discussed by [|^] and who have shown that once the cor- 
rection for this bias is applied, the real intrinsic fraction of Compton-thick AGN is around 20-24%. 
In particular, in our previous work ([|0|]) we have adopted a redshift cut (z=0.015 or 60Mpc) in 
our complete sample of INTEGRAL AGN in order to remove the bias and to probe, although only 
locally, the entire AGN population. Following this reasoning and using the entire AGN sample, 
we should be able to expand our previous analysis and to confirm our initial hypothesis, having 
in mind that the present sample although enlarged is not complete. To allow a comparison with 



the work of [20], we have divided the sample in the same redshift bins (up to z=0.57 considering 
only those AGN with Log L2o-iookeV <46) and plotted the fraction of absorbed (Nh >10^^ cm^^) 
objects compared to the total number of AGN in these bins. The result is shown in Figure ^ (left 
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panel), for the enlarged (red points) and complete sample (black points). A number of consid- 
erations can be made from this figure. First, the bias in z is still present as we keep observing 
a trend of decreasing fraction of absorbed objects as the redshift increases. Second, in the first 
bin the fraction of absorbed objects remains the same as found before: in particular we find that, 
out of 66 objects, 53 (or 80%) have a column density >10^^ and 11 (or \1%±3%) are Compton 
thick. Taking into consideration the fact that the entire sample is not complete, these results are in 
close agreement with those found previously and confirm the original suggestion that our survey is 
able to pick up all AGN, even the most absorbed ones, but only in the local Universe. Finally, the 
fraction of absorbed (and Compton thick) AGN has increased in the second bin from 35% to 57% 
(including 4 Compton thick sources), implying that as the INTEGRAL survey progresses, we are 
able to pick up more absorbed objects among those which are distant and faint and therefore lost 
in previous catalogues. 

We thus conclude that the bias affecting deep hard X-ray surveys of AGN is real but negligi- 
ble if we deal with objects located in the nearby Universe, where the "true" fraction of absorbed 
and Compton thick objects can be estimated with some precision at 80% and ~17% respectively. 
Similar indications are now emerging also from X-ray observations of AGN samples selected in 
other wavebands. For example [||] studied the AGN present in the extended IRAS 12 micron 



sample | |28| ] for which an XMM-Newton observation is available and found that the fractions of 
obscured and Compton thick objects are 62 ±5% and 20±4% respectively. [|l| presented instead 
the XMM-Newton spectral analysis of all 38 Seyfert galaxies in the Palomar spectroscopic sample 



of [15]: they estimate the fraction of absorbed nuclei to be 75±19%, while that of Compton-thick 
sources to be 15-20%. So the numbers obtained from samples selected in various wavebands seem 
to converge at least for what concern the local Universe, while at high z the situation is by far less 
certain. 

4. Hard X-ray spectra 

One problem effecting the study of the hard X-ray spectra of AGN is the low statistical signifi- 
cance of the data which often limits the measurements of parameters such as the reflection fraction 
and the high energy cut-off. A way to overcome this problem is to combine data from INTE- 
GRAL/IBIS and Swift/BAT. Before doing so, it is important to cross calibrate the two instruments: 



one way is the classical method of using the Crab [ ]16| ] , another is to use average spectra of large 
samples of objects belonging to the same class. About 80 objects of the 87 listed in the complete 
INTEGRAL sample have also BAT data available from the 58-month survey ^. A simple power 
law fit to the combined IBIS/BAT points of all these AGN provides a cross-calibration constant of 
I.I4ib0.04, sufficiently close to one to justify the joint use of both datasets. The simple power law 
turns out to be a good description of the 20-100 keV spectra in most cases, with only 35% of the 
sources requiring a more complex fit. The average photon index is around 2 for the entire sam- 
ple and for Seyfert 1 and 2 individually; only Blazars display a slightly flatter spectrum (r'--'1.7). 
Interestingly this average spectrum (especially considering Seyfert 1 alone) is steeper than that 
measured in the CAIXA sample of bright type 1 AGN where r=1.73±0.09 when reflection is in- 
cluded in the fit [Q] ; this could be taken as evidence for the presence of a high energy cut-off in 



http:http://heas arc. gsfc.nasa.gov/docs/swift/results/bs58mon/ 
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Figure 2: Left Panel: Fraction of absorbed objects compared to the total number of AGN as a function of 
redshift from enlarged (red points, and complete (black points, [^o[]) samples respectively. Right Panel: 
pexrav model applied to the combined IBIS/BAT data of 3C 273 with the high energy cut-off fixed to 
IMeV. 



the IBIS/BAT waveband. The presence of this (and the reflection) component can be tested via the 
pexrav model applied to those sources where a simple power law is not a good fit to the data; 
although the improvement is significant in most cases, good constraints on all spectral parameters 
are obtained only for NGC 4151 and IC 4329A. This indicates that, if one wants to get some spec- 
tral information from the hard X-ray spectra alone, some assumptions must be made on at least one 
of the 3 parameters involved in the spectral analysis (photon index, cut-off energy and reflection). 
For example in IGR J21247-I-5058, if we fix the reflection to zero as found by[22, 32], the photon 
index and the cut-off energy can be well constrained (r=1.7^Q } and Ec=155l3y keV) and provide 
values compatible with previous studies. On the other hand, if in the Circinus Galaxy the high 
energy cut-off is set to be around 50keV as measured in the past (e.g. [34]), the photon index and 
the reflection can be estimated with some precision (r=1.80^Q Q3 and R=0.591q 24). Finally in the 
case of 3C 273, where the cut-off is known to be around 1 MeV, the reflection component can be 
constrained to be below 0.3 (see Figure ^ right panel). Unfortunately, such assumptions are not 
possible in all cases and so one has to rely on broad-band data to obtain more stringent constraints 
on source parameters. 



5. Broad-band analysis 

A further step towards the comprehension of hard X-ray selected AGN properties is provided 
by broad-band X-ray studies. In a series of papers [23, 24, 25, 9|], we have combined IBIS spectra 



with data obtained at soft X-ray energies by various satellites such as XMM-Newton, Swift, Chan- 
dra etc., to obtain information on various spectral features for most AGN in the complete sample. 
Here, we report a compendium of the results obtained in these papers. A first result is the great 
variety in shapes observed and the extreme complexity found when modelling most sources; many 
spectral components have to be considered, although they are not present all the time with similar 
strength. The soft excess for example is observed in a good fraction of Seyfert 1 , but dominates 
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dominates only in few of these AGN; it is always present in Seyfert 2, but in 6 objects has a more 
complex shape than a simple scattered power law. Absorption is measured in all type 2 objects, 
while complex absorption is found in 25% of type 1 sources: in this last case one or more layers 
of absorbing material partially covering the source are often found. The average photon index is 
typically 1.7, flatter than the canonical value measured in the hard X-ray band but compatible with 
that measured in the CAIXA sample. Figure || shows the photon index distribution obtained for 
type 1 and type 2 AGN separately: the two distributions as well as the average Y (1.74, a=0.2 and 
1.68, a=0.3 for Seyfert 1 and 2 respectively) are very similar, which immediately suggests that the 
production mechanism is the same in both types of AGN. 

Although high energy spectral components (reflection and the high energy cut-off) are not easy 
to measure due to the sources spectral complexity and non simultaneity of the soft and hard X-ray 
data, they have nevertheless being measured in a number of objects. As is evident in figure || (right 
panel), the distribution of measured cut off energies clusters around 100 keV while the bulk of the 
lower limits on this parameter are found below 300 keV. One must then conclude that high energy 
cut-off is present and cannot be ignored (as often happens) in broad band fits. This result is in line 
with a number of previous studies [|^, ||, ^ but is only barely consistent with those obtained by 
[27], who locate the cut-off energy above 200-300 keV. We also note that our result is in line with 
the synthesis models of the cosmic diffuse background (CXB) which often assume an upper limit 
of ~200 keV for the AGN mean cut-off energy. This choice is basically driven by the intensity and 
shape of the CXB spectrum above the peak, which cannot be exceeded; even a value of 300 keV 
has difficulties in accommodating all available observations and CXB measurements [p^]. 

Finally, Figure ^ (left and right panels) shows a compendium of the results obtained for the 
reflection component in both type 1 and type 2 AGN; note that R is linked to the solid angle D. 
as R=n/27r and that measured values and upper limits are shown separately in the figure. The 
distributions of R are similar in both types as well as the average reflection values; it is also evident 
that large values of R are found in highly variable sources and simply reflect the limitation due to 
the use of non simultaneous data. 



6. Future prospects and conclusions 

The number of AGN selected in the hard X-ray band is continuously increasing thanks to both 
the INTEGRAL and Swift surveys and many of the initial expectations have been fulfilled but also 
many questions have been raised by the observations made so far. Probably the shift in the next 
few years will be to multi-waveband studies of these objects which provide an alternative view due 
to their hard X-ray selection: overall IBIS/BAT objects are less affected by absorption than sources 
selected in other wavebands and are representative of the population of high luminosity AGN. 
Thus multiwaveband studies of them can provide new clues on how AGN work. For example, 
by combining radio and hard X-ray information on a well defined sample of AGN, one can study 
the jet-disk connection in supermassive black holes (Panessa et al. these proceedings) or test if 



jets are made by e+-e^ pairs [13]. Also by surveying hard X-ray selected AGN for water maser 
emission, one can expect to probe the physics of accretion disks, especially in those AGN which 
are often obscured at optical/UV and even X-ray wavelengths by large column densities of gas and 
dust along the Une of sight [^. Already the discovery of a water maser feature in the nucleus of 
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Figure 3: Left Panel: Distribution of photon indices of Seyfert 1 (including Narrow Line Seyfertl) at the 
top and Seyfert 2 at the bottom, as obtained from broad band X-ray data. Right Panel: Distribution of cut 
off energies (measured values at the top and lower limits at the bottom) for the complete sample of AGN as 
obtained from broad band X-ray data. 
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Figure 4: Left Panel: Distribution of the reflection component R for the Seyfert 1 in the complete sample 
as obtained from broad band X-ray data; upper limits on R are shown as a diagonal histogram. Right Panel: 
Same as in left panel but for Seyfert 2. 



IGR J16385-2057 has opened a new field of maser studies in Nairow Lines Seyfert 1 galaxies and 
highlighted the possibility of detecting this emission in elliptical galaxies As underlined 

by [11], modelling of the infrared spectra of a complete sample of hard X-ray selected AGN can 
provide the only way to measure the intrinsic distribution of torus covering factors and thus provide 
the ultimate test for the clumpy torus model. Finally, one should note that BAT/IBIS samples of 
AGN show a large fraction (20-30%, [p7||) of systems in interaction, mergers, dual objects etc. 
Combined radio/optical/X-ray observations of these systems offer a unique opportunity to study 
AGN activation (see for example [18]), probe the relation between interaction and absorption and 
ultimately understand how black holes grow. 

This is of course only the beginning of hard X-ray studies of AGN, since INTEGRAL is 
approved up to 2014 and hopefully beyond, Swift/BAT is still in orbit and performing well and 
Nustar has just joined the group with great expectations, but as Plato would say "The beginning is 
the most important part of the work". 



7 



AGN at hard X-ray energies 



L. Bassani 



References 

[ 1 1 Akylas A. & Georgantopoulos I., 2009, A&A, 500, 999 
Beckmann V., Soldi S., Ricci C. et al. 2009, A&A, 505, 417 
Bianchi S., Guainazzi M., Matt G. et al., 2009, A&A, 495, 421 
Bird A.J., Bazzano A., Bassani L. et al., 2010, ApJS, 186, 1 
BrightmanM. & Nandra K., 2011, MNRAS, 413, 1206 
Burlon D., Ajello M., Greiner J. et al., 201 1, ApJ, 728, 58 
Castangia P., Tarchi A., Panessa F. et al., 201 1, Extremesky 201 1 
Dadina M. 2008, A&A, 485, 417 

De Rosa A., Panessa R, Bassani L. et al., 2012, MNRAS, 420, 2087 
Elitzur M. & Shlosman I., 2006, ApJL, 648, 101 
Elitzur M., 2012, ApJL, 747, 33 
Filippenko A. V. & Ho L. C, 2003, ApJL, 588, 13 
Ghisellini G. 2012, MNRAS, 424, 26 
GilU R., Comastri A. & Hasinger G. 2007, A&A, 463, 79 
Ho L. C, Filippenko A. V., Sargent W. L. et al., 1997, AAS, 189, 122.09 
Kirsch M. G. et al., 2005, SPIE, 5898, 22 

Koss M., Mushotzky R., Veilleux S., Winter L., 2010, ApJ 746, 22 
Koss M., Mushotzky R., Treister E. et al., 2012, ApJ, 716, 125 
Krivonos R., Revnivtsev M., Lutovinov A. et al., 2007, A&A, 475, 775 
Malizia A., Stephen J. B., Bassani L. et al., 2009, MNRAS, 399, 944 
Malizia A., Bassani L., Bazzano A. et al. 2012, MNRAS, 426, 1750 
MoUna M., Giroletti M., Malizia A. et al., 2007, MNRAS, 382, 937 
Molina M., Bassani L., Malizia A. et al., 2009, MNRAS, 399, 1293 
Panessa P., Bassani L., de Rosa A. et al. 2008, A&A, 483, 151 
Panessa P., de Rosa A., Bassani L. et al. 201 1, MNRAS, 417, 2426 
Perola G. C, Matt G., Cappi M. et al., 2002, A&A, 389, 802 
Ricci C, Walter R., Courvoisier T. J.-L. et al., 201 1, A&A, 532, 102 
Rush B., Malkan M. A. Spinogho L., 1993, ApJS, 89, 1 
Silverman J. D., Green P J., Barkhouse W. A. et al., 2005, ApJ, 618, 123 
Tarchi A., Castangia P, Henkel C. et al., 201 1, A&A, 525, 91 
Tarchi A., Castangia P, Columbano A. et al., 201 1, A&A, 532, 125 
Tazaki P, Ueda Y., Ishino Y. et al., 2010, ApJ, 721, 1340 
Urry C. M. & Padovani P, 1995, PASP, 107, 803 
Yang Y, Wilson A. S., Matt G. et al., 2009, ApJ, 691, 131 
Winkler C, Courvoisier T.J.L., Di Cocco G., et al., 2003,A&A, 41 1, LI 



8 



